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We propose and experimentally demonstrate detecting small single-cycle and few-cycle signals by 
using the symmetric double-well potential of a radio frequency superconducting quantum interfer¬ 
ence device (rf-SQUID). We show that the response of this bistable system to single- and few-cycle 
signals has a non-monotonic dependence on the noise strength. The response, measured by the 
probability of transition from initial potential well to the opposite one, becomes maximum when the 
noise-induced transition rate between the two stable states of the rf-SQUID is comparable to the sig¬ 
nal frequency. Comparison to numerical simulations shows that the phenomenon is a manifestation 
of stochastic resonance. 


It is a long-held belief that noise is detrimental or even 
destructive to detecting signals which often appear as 
weak periodic modulations. However, during the last 
35 years theoretical and experimental investigations have 
shown that in nonlinear systems a proper amount of noise 
can actually increase the signal-to-noise ratio (SNR) and 
thus become beneficial for signal detection. This interest- 
ingphenomenon is named as stochastic resonance (SR) 
[2-1^. For example, suppose that a particle is moving in 
a periodically perturbed symmetric double-well poten¬ 
tial under the influence of a Gaussian-white noise such 
as thermal fluctuation. Then SNR of the power spec¬ 
tral density of the particle’s trajectory displays a broad 
maximum when the rate of inter-well transitions, which 
depends on noise strength exponentially, is comparable 
to the frequency of periodic signal. This is the essence of 
SR. 

Due to its simplicity and ubiquity of the underlying 
mechanism, SR has attracted much interest from physi- 
cists, chemists, biologists, and electronic engineers |ll-[lll|. 
It has also been observed in Josephson junction-based 
systems (T2l - [l^ . which have recently attracted much in¬ 
terest and been app lied in many fields such as quantum 
information [l7H19|. However SR has been only investi¬ 
gated for periodic signals that last many cycles. Namely, 
only the steady-state properties of the noisy periodically 
driven systems have been studied. On the other hand, 
in a variety of science and engineering disciplines, it is 
a significant challenge to detect small signals which not 
only last a few cycles but also are buried in noise. Up 
to this point, whether SR can also enhance single-cycle 
signal detection remains an open question. 

In this Letter, we report on the observation of SR in 


*Electronic address: gzsun@nju.edu.cn 


a radio frequency superconducting quantum interference 
device’s (rf-SQUID’s) [lO, [HI response to weak single¬ 
cycle and few-cycle signals by measuring the inter-well 
transition probability as a function of the noise strength 
D and the signal frequency fs systematically. Our ex¬ 
perimental and numerical results show that one can dis¬ 
till small single-cycle and few-cycle signals from noisy 
environment by using bistable systems configured as bi¬ 
nary threshold detectors. The maximum sensitivity is 
achieved at the value of D that matches well with the 
position of SR. We also show that the sensitivity of de¬ 
tecting single-cycle signals is similar to that of detecting 
many-cycle signals. 

In our experiment we use an rf-SQUID, which is a 
superconducting loop of inductance L interrupted by a 
Josephson junction of critical current Ic, as our bistable 
detector. An optical micrograph of the sample is shown 
in the inset of Fig. 1. The Josephson junction is made of 
Nb/AlOx/Nb on a silicon substrate. The critical current 
Ic and the capacitance C of the junction are approxi¬ 
mately 0.80 ^A and 90 fF, respectively. The inductance 
L of the Nb superconducting loop is approximately 1053 
pH. The potential energy of the rf-SQUID is given by 
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where d)o is the flux quantum and Ej = /cd)o/27r is the 
Josephson coupling energy of the junction. The shape of 
the double well potential can be controlled in situ by a 
flux bias d>e applied via a flux bias line coupled induc¬ 
tively to the rf-SQUID. In particular, at = ‘I’o/2 the 
SQUID has a symmetric double-well potential separated 
by a barrier AUq as shown in Fig. 1. For the SQUID 
studied here, AUo/kB — 12.3 K, where ks is the Boltz¬ 
mann’s constant. The dynamics of the rf-SQUID, identi¬ 
cal to that of a fictitious flux particle of mass C moving 
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in the potential [/(<&) with a friction coefficient R is 
governed by the corresponding Langevin equation: 
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Here, In is the noise current and R is the damping re¬ 
sistance of the Josephson junction. Without externally 
injected noise, /„ and R are related by the fluctuation- 
dissipation theorem )) = — t ) in ther¬ 

mal equilibrium, where T is temperature. The small os¬ 
cillation frequency of the system around the bottom of 
the potential wells is denoted as wo- At T ^ Huio/kB, 
where h is the Planck constant, thermal activation causes 
inter-well hopping with the characteristic transition rate 
given by 
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where a* is a damping dependent constant of order of 
unity. When transitions are dominated by an external 
noise source of strength D ^ fc^T, the denominator in 
the exponent of Eq. m is replaced by D which is pro¬ 
portional to the square of the rms noise current, ^ms) 
applied to the system. For the sake of simplicity, here¬ 
after we set fcs = 1 so that D is measured in units of 
kelvin. Note that because the potential is symmetric, Pq 
is identical for left-to-right and right-to-left transitions. 

Because all key parameters of the rf-SQUID potential 
and its control circuit can be accurately determined, the 
double-well potential of the rf-SQUID is an ideal system 
for investigating SR [l^ - [l^ and noise-enhanced detection 
of single-cycle and few-cycle signals. In our experiment, 
the Gaussian-white noise has a bandwidth of about 9 
MHz, which is generated by an arbitrary waveform gen¬ 
erator. The signal and noise are applied to the rf-SQUID 
through a second flux bias line with higher bandwidth 
(up to 18 GHz). The relationship between D and In,ms 
is calibrated by measuring Pq versus In rms s-iid compar¬ 
ing the result to Eq. 

As shown schematically in Fig. 1, each measurement 
cycle begins by ramping up the quasi-static flux bias from 
0 to $ 0 /2 to prepare the flux particle in the left side of the 
symmetric double-well potential. This is followed by ap¬ 
plying a single-cycle modulation of flux bias that causes 
the potential barrier to oscillate as AC/q -I- Eq sin(27r/st), 
where Eq is proportional to the amplitude of the flux 
modulation which is kept at Eq = O.OTAUo — 0.86 K 
in the experiment. The position of the flux particle is 
measured by using a dc-SQUID switching magnetome¬ 
ter inductively coupled to the rf-SQUID, either after a 
single signal cycle or a fixed duration of signal time as 
discussed later, as a function of 0.5 K< D < 2.0 K and 
10 kHz < fs < 200 kHz. The quasi-static flux bias is 
then ramped down to zero to complete the measurement 
cycle. To obtain the fractional population in the right 
well pb, the procedure is repeated 2000 times at each 
value of D and fs- All data are measured at T Ri 20 mK 


D in a cryogen-free dilution fridge carefully shielded 
from the environmental electromagnetic interference so 
that the effects of thermal fluctuation and extra noise on 
the experiment are negligible. 

We first measure pa as a function of the noise strength 
D by using single-cycle signals e±{t) = ±Eo sin(27r/st) as 
depicted in Fig. 1, where the signal frequency /s = 10 
kHz. The result is shown in Fig. 2(a). The noise strength 
D is varied between 0.5 K and 2.0 K. Therefore, transi¬ 
tions between the two potential wells at Eq = 0 (no sig¬ 
nal) are noise activated. Note that without the noise, 
e{t) alone would be too small to cause transitions be¬ 
tween the potential wells because (AC/q — £o)/T > 500. 
Hence, D > 0.5 K is required for the flux particle to hop 
from one well to the other within the duration of each 
signal cycle. On the other hand, at Eq = 0 the transi¬ 
tion rate Tp grows exponentially from approximately 1 /s 
at D = 0.5 K to greater than 10^/s at D = 2.0 K, as 
shown in the inset of Fig. 2(a). It can be seen that for 
D < 0.7 K the population of the right well pa is negli¬ 
gible at t = r = l//s. The data indicated by the blue 
squares in Fig. 2(a) are taken with £+{t) which show 
that as D is increased from 0.7 K, pn rises rapidly to 
reach maximum when ro(Dm) 2fs, where D^i denotes 
the noise strength corresponding to the maximum pn- 
When D > D^, the probability of hopping back from 
the right well to the left well increases rapidly, causing 
Pr{D) to decrease. Finally, when D » Eq the population 
of each potential well is equalized to 50%. 

SR has two most prominent signatures: One is the 
peak in the system’s response versus noise strength D. 
The other is the position of the peak and the sig¬ 
nal frequency satisfying ro(Dm) ^ 2fs, or equivalently, 

l/Dn, -ln(/,) according to Eq. ©. In Fig. 2(b), 

where £+(t) is applied, we plot pR versus 1/D and fs 
which shows clearly both signatures of SR. In particu¬ 
lar, the nearly linear relationship between ro(Di„) and 
fs is demonstrated as shown in the inset of Fig. 2(b). 
The slope of ro(Di„) versus fs obtained from the best-fit 
to a line is 2.7, which is consistent with the numerical 
result previously obtained for SR under continuous mod¬ 
ulation (^ . In addition, we numerically calculate the 
power spectral density S{fs) of the flux particle’s trajec¬ 
tories <I>(t) generated by Monte Carlo simulation of Eq. 
([ 2 ]). It is found that S{fs) reaches its maximum at the 
same value of D as pR does. We thus conclude that SR 
plays a central role in the bistable system’s response to 
single-cycle signals. 

In order to compare the result of our measurements 
with that of numerical study over the entire parameter 
space covered by the experiment, we adopt the two-state 
model Bm and introduce the rate equation: 


’^P^ = -T.pR{t)+V+[l-pR{t)] (4) 

with the initial condition pl{Q) = 1, Pr{0) = 0. Here, pr 
{pl = 1 — Pr) is the fractional population of the right 
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(left) potential well. When Eq ^ 0, the barrier height is 
oscillating between AC/q i s-nd the transition rates are 
time-dependent 


r±(t) = To exp 


±eo sin(27r/t) 
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where r+ and r_ denote the rates of left-to-right and 
right-to-left transitions, respectively. Fg is given by Eq. 
©. Using the system parameters given above, we nu¬ 
merically integrate Eq.Q to obtain pR{t) as a function 
of fs and D. The result is shown in Eig. 2(c). It can be 
seen that the key features of the experimental data are 
well reproduced. 

Next, we show that the sensitivity of detecting single¬ 
cycle signals is comparable to that of detecting many- 
cycle signals and that one can predict the population dis¬ 
tribution of the bistable systems at the end of A-cycle 
modulations pr^n = Pr(,Nt) = 1 — pl,n from that of 
single-cycle modulation pr^i. It is straightforward to ob¬ 
tain the recursion relation 


PR,n+l = PL,nP+ + PR,n{^ — P-) ( 6 ) 

= (1 - PR,n)P+ + PR,nO- - P-) ■ 

The first (second) term of the r.h.s. of Eq. (O is the 
fractional population of the left (right) well at t = nr 
that ends (remains) in the right well at t = {n + 1)t. 
P+{-) is the probability of switching from the left (right) 
to the right (left) well during the time interval nr < t < 
(n + 1)t. Notice that with the single-cycle perturbation 
£rj_(0 < t < r) and the initial condition pr o = 0, one 
has P± = pr^i by taking into consideration the spatial 
and temporal symmetry properties of the rf-SQUID po¬ 
tential and e±. Thus, we can obtain P± directly from the 
data presented in Fig. 2(a). Because Eq. ([6]) is valid for 
arbitrary noise strength D and signal frequency fs, we 
can compute pr,n from P± for any integer N > 1. We 
find that as N increases pr^n converges rapidly. In order 
to investigate the dependence of SR on the number of 
signal cycles N, we modify the experimental procedure 
by changing the duration of the applied signal and noise 
from T to 0.3 ms. Thus, we have A = 3 for /g = 10 


kHz, which increases ultimately to A = 60 for fg = 200 
kHz. In Fig. 3(a), the measured pr,n is plotted against 
1/D and fg, which compares well with pr^n computed 
from Eq. (jS]) by using the measured P± as inputs [see Fig. 
3(b)] and that obtained by solving the corresponding rate 
equation (g]) [see Fig. 3(c)]. The results presented in Fig. 
3 all have two distinctive features: (i) The threshold noise 
strength Do which demarcates the blue region {pR ,N ~ 0) 
and the yellow region depends weakly on the number of 
signal cycles N, and (ii) 1/D^ oc — ln(/g) remains valid 
for the entire range of 3 < iV < 60. As shown in the 
inset of Fig. 3(a), the dependance of ro(I?m) on fg is ap¬ 
proximately linear with a slope of about 2.6. These two 
features strongly indicate that the sensitivity of detect¬ 
ing single-cycle signals is similar to that of many-cycle 
and continuous wave signals and that SR does exist in 
the systems driven by small single-cycle signals. 

In summary, using an rf-SQUID as a prototypical 
bistable system, we have demonstrated the existence 
of SR with single-cycle perturbation to the symmetric 
double-well potential of the system. Furthermore, we 
have investigated the possibility of exploiting SR for de¬ 
tecting small single-cycle and few-cycle signals in noisy 
environment. We have found that a proper amount of 
noise can lead to SR which enhances the sensitivity of 
detection. Our work provides insights into the behav¬ 
ior of bistable systems under the combined influence of 
weak single-cycle (or few-cycle) periodic modulation and 
noise. Because conventional techniques, such as phase 
sensitive lock-in and heterodyne detection schemes, are 
not applicable to detecting single-cycle and few-cycle sig¬ 
nals buried in noise, the method demonstrated here is 
promising for applications where signals are unavoidably 
mixed up with noise and only last a very small number 
of cycles. 
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FIG. 2: (a) Measured pr{D) with a single-cycle sinusoidal 
signal, fs = 10 kHz. Data indicated by the blue squares and 
red circles correspond to e+{t) and e-{t), respectively, which 
agree with the numerical calculation (black lines). The in¬ 
set shows the linear dependance of In(ro) on 1/D. (b) pu 
as a function of 1/D and fa with single-cycle sinusoidal sig¬ 
nals e+{t) = eo sin{2n fat) . As shown in the inset, when pR 
reaches a maximum, ro(Din) ~ 2/s, where Dm denotes the 
noise strength corresponding to maximum pR. These results 
are consistent with the theoretical hypothesis of SR. (c) Nu¬ 
merical calculation of pR as a function of 1/D and fa with 
single-cycle sinusoidal signals e+{t) = eo sin(27r/st), which 
agree well with the experimental data. 


FIG. 1: A time prohle of manipulation and measurement. Po¬ 
tential wells at several key moments are also plotted. The in¬ 
set shows an optical micrograph of a Nb/AlOx/Nb rf-SQUID 
with an inductively coupled dc-SQUID and flux bias lines. 
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FIG. 3: (a) Measured pn as a function of 1/D and with 
a constant duration (0.3 ms). SR remains as shown in the 
inset, (b) pn as a function of 1/D and fa with a constant 
duration (0.3 ms), derived from the recursion relation Eq. (|6]) 
and P± obtained from the experimental results when single¬ 
cycle signals are used, (c) Numerical calculation of pR as a 
function of 1/D and fs with a constant duration (0.3 ms). 









